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Abstract - Glass formers exhibit, upon an oscillatory excitation, a response function whose 
imaginary and real parts are known as the loss and storage moduli respectively. The loss modulus 
typically peaks at a frequency known as the a frequency which is associated with the main 
relaxation mechanism of the super-cooled liquid. In addition, the loss modulus is decorated by 
a smaller peak, shoulder or wing which is referred to as the /3-peak. The physical origin of this 
secondary peak had been debated for decades, with proposed mechanisms ranging from highly 
localized relaxations to entirely cooperative ones. Using numerical simulations we bring an end 
to the debate, exposing a clear and unique cooperative mechanism for the said /3-peak which is 
distinct from that of the «-peak. 



Introduction: In Fig. Q] one observes a typical loss 
1 modulus, in this case of a film of the metallic glass 
Zr65Al7.5Cu27.5 which was forced at a fixed frequency 
oj = 5440Hz while the temperature was changed from 
450K to 850K pQ. The main peak represents the typical a 
relaxation which is the slowest mode of relaxation which is 
typical to glass formers. When the temperature was at a 
value such that this relaxation frequency was at resonance 
with the fixed forcing frequency, the loss of energy regis- 
tered a maximum. Upon reducing the temperature one 
sees a wing which becomes obvious at about 650K, this 
wing is associated with a typical frequency that is higher 
than the a frequency. This is referred to as a secondary 
relaxation or the (3 wing, and in different glassy materials 
it appears more or less prominent, sometime as a separate 
peak, and sometime shoulder or a wing [2]. Note at 
this point that the literature is full of other phenomena in 
glassy relaxation that are referred to as f3 [i 8 . In this 
Letter we only deal with the one presented here as in Fig. 
[TJ This is also referred to as the Johari- Goldstein process 
or the slow /3-process [3]. 

The existence of the /3-peak was known for a long time, 
but its physical origin remained debated. Past explana- 
tions can be roughly grouped into two main groups, those 
involving an extension of the a-relaxation process and 
those proposing an independent relaxation mechanism. In 
the first group one can quote for example Dixon et al. [9] , 



10 



Zr Al Cu ■ film 
(i relaxation (HN-function) / f 

t 

wing / 




450 500 550 600 650 700 750 800 850 
T[K] 



Fig. 1: Color online; a typical loss modulus as a function of 
temperature. Shown are experimental results for a film of the 
metallic glass Zr65 AI7.5CU27.5 which was forced at a fixed fre- 
quency uj — 5440Hz. Note the well defined a peak and the 
shoulder, or wing, which appears at lower temperatures and is 
referred to as the j3 peak. [TJ.Havriliak-Negami (HN) curve is 
plotted using parameter a = 1 and 7 = 0.6. Please see Ref. [I] 
for the form of the HN function and other details. 
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whereas the second group is exemplified by |10j . There 
seems to be no consensus even on the question whether 
the /3 process is local or cooperative. The aim of this Let- 
ter is to disperse the fog that had accumulated over the 
physics of the /3-peak and to reach a unique interpreta- 
tion of the relevant physics. To this aim we use numerical 
simulations as described next. Of course, numerical simu- 
lations are limited in their ability to probe long times, and 
we arc therefore constrained to examine temperatures for 
which the increase in glassy relaxation times is only 3 to 

4 orders of magnitude compared to twelve to fourteen or- 
ders in experiments. Nevertheless, we do not expect a 
fundamental change in physics as the system cools down 
towards longer and longer relaxation times. This is the 
basic assumption on which our conclusions rely. 

Numerical simulations: As a first step we have re- 
produced the phenomenology of the /3-peak on the ma- 
chine. To this aim we employed as an example one of the 
standard models of glass formers, i.e. a binary mixture 
of point particles interacting via Lennard- Jones poten- 
tial with three different characteristic interaction lengths 
a ss = l,cr« = 1.4 and a 8 £ = 1.18. Details of the poten- 
tials can be found for example in Ref. pj] . Of course the 
phenomenon discussed appears in other models in much 
the same way, as we have checked by running similar sim- 
ulations on other models. In the present simulations the 
unit of length is a ss , temperature is measured in units 
such that the Boltzmann's constant equals unity, and time 
in units of inter-particle distance over the mean velocity. 
Two dimensional samples, each containing 4900 particles, 
were firstly quenched slowly (reducing the temperature by 

5 x 10~ 4 in every step) to an inherent state configuration 
at T = using a gradient energy method. Wc employ 
periodic boundary conditions keeping the volume of the 
box fixed. In a second step the samples were raised to 
a chosen temperature T and were subjected to an oscil- 
latory shear strain 7 = 70 sin(wt) according to the affine 
transformation 

Xi ->• Xi + {jdtyyi, yi ->■ yi . (1) 

in addition we impose Lees-Edwards boundary conditions. 
The maximum affine strain is 70 = 0.005, and Xi, yi are the 
position of the i'th particle. The Gear predictor-corrector 
algorithm was used to integrate the shod equations of mo- 
tion p] . 

The stress in the system was measured using 50 different 
realizations, and the mean shear stress was fitted to an 
oscillatory function r xy — tq sm(u>t+S). The S corresponds 
to the phase shift between the strain and the stress. Now, 
we can calculate the dynamic shear modulus, 

G = —e lS (2) 
7o 

where G = G + iG is a complex number that depends 
on the temperature and the shearing frequency. G and 
G arc known as the storage and the loss moduli. 
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Fig. 2: Color online; the loss modulus G as a function of 
temperature from numerical simulations (w = 0.01). Observe 
the clearly defined a peak and the distinct ft wing. 



To produce data which are analogous to the experimen- 
tal ones in Fig. [T] we studied G when the system was 
subject to such a simple oscillatory shear with a fixed fre- 
quency fi for a range of temperatures. A typical result 
is shown in Fig. [2] which pertains to uj = 0.01 and T in 
the range 0.5< T <6. In a satisfactory agreement with 
the experimental picture we observe a clear a-peak deco- 
rated with a distinct /3-wing in the range of temperatures 
1.0< T <2.0. Similar simulations done at a fixed tempera- 
ture with varying frequencies resulted in the same evidence 
for the existence of the /3-process as seen in experiments. 
We thus feel encouraged to proceed to study the nature of 
this process. 

The /3-process, local or cooperative?: presently we 
can use the power of numerical simulations to answer the 
first crucial question, is the /3-process local or cooperative? 
To this aim we have selected randomly a small fraction 
of the particles in our system, and pinned them to the 
coordinate frame that moves with the affine transforma- 
tion, not allowing them to participate in the relaxation 
dynamics. We expect that such a procedure should cer- 
tainly affect the a-peak which is known to be cooperative 
[I21H3]. As for the /3-wing, if it is associated with a lo- 
cal process, a small fraction of pinned particles should not 
make a big difference. It turns out that a small fraction 
of pinned particles, say larger than 2.5%, is sufficient to 
suppress the /3-wing altogether. In Fig. [3] we show the 
results of this exercise with varying percentage of pinned 
particles, from 0% to 10%. While the a-peak is affected as 
expected, slightly suppressed and pushed to higher tem- 
peratures (the decay process is pushed to lower frequen- 
cies) |6], the /3-wing is suppressed altogether when the 
concentration of the pinned particles exceeds 2.5%. We 
can thus conclude that the /3-process must be coopera- 
tive rather than local. This immediately excludes the first 
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Fig. 3: Color online; Upper panel: the loss modulus as shown in 
the previous figure, and the loss modulus when a small fraction 
of particles get pinned. The consequence of varying concentra- 
tions of pinned particles are shown. Note that it is sufficient to 
pin more than 2.5% of the particles to eliminate the /3 wing al- 
together. Lower panel: collapse of the data shown in the upper 
panel, obtained by rescaling the temperature with a character- 
istic G a and T a , the value and the position of the a peak which 
depends on the concentration of pinned particles. We see that 
the /3 peak is eliminated with all the concentrations larger than 
2.5%. The grey curve is a Gaussian fit to the remaining a pro- 
cess. 



group of mechanisms discuss above. 

Which cooperative mechanism then?: to un- 

fathom the nature of the cooperative (3 process, we exam- 
ined carefully the dynamics in our system on time scales of 
the order of the cage dynamics, i.e. for time scales of the 
order of unity in our re-scaled units. In Fig. [3] we show the 
self part of the intermediate scattering function F s (q, t) for 
a fixed temperature and various pinning density. F s (q,t) 
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Fig. 4: Color online. The self part of the intermediate scat- 
tering function for T = 3 for several pinning densities in a 2D 
system. Note the change of the shape of the function in the 
encircled area. 



is defined as 

P '(9'*)4(S e,fW(,K,{0)1 )' (3) 

where N' is the number of particles which are mobile, 
(. . .) represent thermal average and rj(t) is the position 
vector of particle j at time t. q is chosen uniformly on 
the surface of a sphere of radius q = \q\ = 27r/ro, where 
r = 2 1 / 6 is the position of the minimum in the potential 
function as a function of distance. This value roughly cor- 
responds to the first peak position in the static structure 
factor. We took 128 different directions for the averaging 
of F s (q, t) . We observe a strong slowing down in the dy- 
namics as we increase the pinning density: the effective r a 
becomes larger. But we also observe a hump in the curve, 
on time-scales corresponding to the (3 process, which al- 
most disappear when the concentration of the pinned par- 
ticles vanishes. This hump can be explained by a recoil 
motion of a free particle when it collides with a pinned par- 
ticle. Evidently this phenomenon is less pronounced when 
the collision occurs with two free particles. In this case, 
one particle can transfer momentum to its neighbor etc, a 
process that evolves into a collective motion of number of 
particles due to the collisions. This cooperative motion is 
suppressed by the pinned particles. 

In left panel of Fig. [5] we show the cumulative displace- 
ment of all the particles after four time units (i.e. at t = 4 
with respect to their position at t = 0). We observe long 
chains of correlated motion of particles all occurring at 
thermal equilibrium which is maintained only by the in- 
fluence of random forces. We emphasize that no external 
forces were applied. The darker arrows are to indicate the 
particles that move more than 40% of the typical inter- 
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Fig. 5: Color online: Left Panel: a graphic representation of the cooperative motion that is associated with the f3 wing. Note the 
chains of particles that have moved coherently during a time span of four time units. The particles that move more than 40% 
of the typical inter-particle distance are marked in dark blue. Right Panel: similar graphic representation of the suppression of 
the majority of the cooperative motion that is responsible for the j3 wing by the addition of 2.5% pinned particles. In contrast 
to the previous figure, here one needs to look at cumulated motions for 15 time units to see the remnant correlated motion. 



particle distance. In comparison, with 2.5% of pinned 
particles, there is hardly any such correlated motions on 
the scale of 4 units of time. One has to increase the time 
window to about 15 units of time, cf. the right panel of 
Fig. [5l in order to observe the remnant weak cooperative 
motion. 

In summary, we have shown here that the (3 wing mea- 
sured in loss moduli has to do with a cooperative motion 
which is however associated with very fast time scales, of 
the order of the cage time. The cooperative motion con- 
sists of chains of particles that move cooperatively when 
the system is forced in an oscillatory fashion. Thus co- 
operative yes, but an extension of the a process no. We 
trust that the clear cut numerical experiments discussed 
above should remove any doubt as to the origin of the (3 
wing in the loss modulus. 
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